Quantum Phases of Cold Polar Molecules in 2D Optical Lattices 
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We study the quantum phases of hard-core bosonic polar molecules on a two-dimensional square 
lattice interacting via repulsive dipole-dipole interactions. In the limit of small tunneling, we find 
evidence for a devil's staircase, where Mott-solids appear at rational fillings of the lattice. For finite 
tunneling, we establish the existence of extended regions of parameters where the groundstate is a 
supersolid, obtained by doping the solids either with particles or vacancies. We discuss effects of 
finite temperature and finite-size confining potentials as relevant to experiments. 
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Trapped atomic and molecular quantum gases al- 
low the realization of quantum lattice models of strongly 
interacting particles The preparation of quantum 
gases of groundstate polar molecules with strong electric 
dipole moments opens the way to the study of lattice 
models with tunable long-range interactions in atomic- 
molecular-optical (AMO) setups [HHf- In this context, 
the challenges are to (i) identify the fundamental Hamil- 
tonians underlying the physics of strongly interacting 
dipolar gases, (ii) analyze the associated quantum phases 
in the context of their realization in AMO setups. This 
entails (iii) studying the preparation and detection of 
these phases in the presence of finite-size confining poten- 
tials and finite temperature, as relevant to experiments. 

In this letter we analyze the microscopic tight- 
binding Hamiltonian realizable with polar molecules 
trapped in optical lattices under collisional stability. This 
is a 2D Hubbard-likc model for hard-core particles with 
infinite-range interactions which we show to display novel 
quantum phenomena with no counterpart in the atomic 
case. These include: (a) a devil's staircase (DS) of Mott 
solids at rational lattice fillings, and (b) supersolid phases 
(SS) obtained by doping solids with either vacancies or 
additional particles. The infinite-range interactions sta- 
bilize the SS as the low-energy phase for a large range of 
system's parameters, raising interesting prospects for its 
realization with polar molecules. While various kinds of 
SS phases (but not the DS) have been found in mod- 
els with shorter-range interactions [e.g., nearest- 
neighbor (NN) interactions with soft-core particles, or 
NNN interactions], here the emphasis is in simulating 
experimentally relevant systems with up to N ~ I0 3 par- 
ticles, at finite temperature, and for finite-sizes with har- 
monic confinement. Results for experimental observables 
with current AMO technology are presented. 

We consider N bosonic polar molecules aligned by a 
static electric field with induced dipole moment d = \f~D, 
implying strong dipole-dipole interactions. Collisional 
stability is reached confining the molecules to a 2D plane, 
using a strong transverse trapping field, e.g a ID opti- 
cal lattice, with a harmonic oscillator frequency lj±_ to 



prevent collapse due to attractive forces between aligned 
dipoles, for intcrparticle distances larger than a m i n = 
(12D/ma; 2 _) 1 / 5 , with m the mass of a molecule Q. An 
additional 2D optical lattice with spacing a > a nlln and 
harmonic oscillator frequency oj confines the particles in- 
plane. The following single-band Hamiltonian for hard- 
core bosons on a 2D square lattice is then obtained pro- 
vided tiuj^ ^> frw > {D/a 3 ,kBT}, with T the temper- 
ature, and with the requirement that the initial system 
has no doubly occupied sites Q 



The first and second terms in Eq. (JXJ) describe the stan- 
dard kinetic energy with hopping rate J and the repul- 
sive dipole-dipole interaction with strength V = D/a 3 
and Tij = \i — j\, respectively; hi and b\ are bosonic 
operators with b\ 2 = and = b\bi, fi is the chem- 
ical potential and ft = moj 2 a 2 /2. For a gas of RbCs 
polar molecules (d = 1.25Debye) with transverse and in- 
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FIG. 1: (color online) (a) Phase diagram of model JT]) with 
Q, = and at T — 0. Lobes: Mott solids (densities indicated); 
SS: supersolid phase; SF: superfluid phase. DS: devil's stair- 
case. Panels (b-d): sketches of the groundstate configuration 
for the Mott solids in panel (a), with density p — 1/2, 1/3 and 
1/4, respectively. 
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plane trapping Vo.±/Er — 40 and Vq/E-r = 4. respec- 
tively, and lattice spacing a = 400nm > a m ; n ~ lOOnm, 
uj ± /2ir > cj/2tt > D/(a 3 h), with D/{a 3 h) ~ 3.5kHz, 
and thus the single-band Hamiltonian Eq. (JTJ) is valid, 
provided T < 200nK. With an in-plane tunneling rate 
J/h ~ 120Hz, the ratio J/V can be tuned to any value 
J/V > 0.03 by varying the strength of the DC field. 
We have studied the quantum phases of Eq. (JTJ by 
means of large scale Monte-Carlo simulations based on 
the Worm Algorithm Q , using no cutoff in the range of 
the dipolc-dipole interaction [lOj. While the focus is on 
experimentally-relevant trapped systems with finite T, 
we find it convenient to first discuss the phase-diagram 
in the homogeneous situation, and then use these results 
to explain the phases for ^ 0, and observables in ex- 
periments. This is followed by more detailed discussions 
of some aspects of the various phases. 

Homogeneous case: Our exact results for = are 
summarized in Fig. [T^a), the computed T — phase- 
diagram as a function of p and J, with f < p/V < 6 and 
J/V > 0.02 (unshaded area). For small-enough hopping 
J/V <C 0.1, we find that the low-energy phase is incom- 
pressible (dp/ dp = 0, with p the filling- factor) for most 
values of p. This parameter region is labeled as DS in 
the figure and corresponds to the classical devil's stair- 
case, i.e. a succession of incompressible ground states, 
dense in the interval < p < 1, with a spatial structure 
commensurate with the lattice for all rational fillings [ll[ , 
with no analogue for shorter-range interactions. For fi- 
nite J, three main solid Mott lobes emerge with p = 1/2, 
1/3, and 1/4, named checkerboard (CB), stripe, and star 
solids, respectively. The corresponding groundstate con- 
figurations are sketched in panels (b-d) . These large Mott 
lobes are found to be robust in the presence of a confining 
potential and finite T (sec below), and thus are relevant 
to experiments. We notice that the shape of the solids 
with p = 1/2 and 1/4 away from the tip of the lobe can 
be shown to be qualitatively captured by mean-field cal- 
culations, while this is not the case for the stripe solid at 
filling 1/3 which has a pointy- like structure characteristic 
of fluctuation-dominated ID configurations. Mott lobes 
at other rational filling factors, e.g. p = 1 and 7/24, are 
present, however not shown here. 

For large enough J/V, the low-energy phase is su- 
pcrfluid, for all p. At intermediate values of J/V , how- 
ever, we find that by doping the Mott solids either with 
vacancies (removing particles) or interstitials (adding ex- 
tra particles) a supersolid phase can be stabilized, with 
coexisting supcrfluid and crystalline orders (we find no 
evidence of SS in the absence of doping). We find that 
the solid/superfluid transition consists of two steps, with 
both transitions of the second-order and the supersolid 
as the intermediate phase (see Fig. |4] below). Remark- 
ably, the long-range interactions stabilize the supersolid 
in a wide range of parameters. For example, a vacancy 
SS is present for densities 0.5 > p > 0.43, roughly inde- 
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FIG. 2: (color online) Spatial density profile in 2D for TV ~ 
1000 particles in a harmonic potential. Phases are indicated 
(CB, SR, ST stand for checkerboard, star, and stripe solids 
respectively), (a-b) V/J = 15, p/J = 55, fi/J = 0.05 
and T/J = 0.0377, with temperature annealing performed in 
panel a); (c) V/J = 5, u/J = 19, fi/J = 0.01 and T/J = 0.1; 
(d) V/J = 20, u/J = 51, fi/J = 0.04 and T/J = 0.25. 



pendent of the interaction strength. For experiments, a 
fundamental question is the observability of the phases 
described above for finite T . In particular, for the SS 
phase we show in Fig. [5] that by increasing T it melts 
into a featureless normal fluid via a two-step transition, 
the intermediate phase being a normal fluid with finite 
density modulations, similar to a liquid crystal. 

Harmonic trap and experimental observables: The 
recent experimental achievement of single-site address- 
ability in optical lattices using electron and optical mi- 
croscopy allows for a direct, in-situ, observation of parti- 
cle positions and particle-particle correlations in experi- 
ments [121 ] . Thus, the key observables for present and fu- 
ture experiments are the in-situ density distribution and 
particle correlations, from which the phases above can 
be detected. The question is how the phases described 
in Fig. [1] will be seen in an experiment. Here, we provide 
snapshots of particle configurations for realistic experi- 
mental situations with N ~ 10 3 particles trapped with 
harmonic confinement, and small finite T. 

In Fig. [2] we show snapshots of the spatial density 
distribution in the lattice (shown is a single quadrant). 
Each circle corresponds to a different site, and its radius 
is proportional to the local density. In panels a) and b) , p 
has been chosen such that particles at the trap center are 
in the CB phase, with very small T. The density profile 
shows a wedding-cake structure, with concentric Mott- 
lobes with density p — 1/2 and 1/4, analogous to the 
shells with contact interactions However, while the 
system parameters are the same in both figures, panel (a) 
shows regular CB and star patterns, while in panel (b) 
extended defects are present in the CB phase and the star 
is barely visible. This is due to the different preparation 
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FIG. 3: (color online) p vs. fi. (a): Solids and SS for a system 
with linear size L — 12 and J/V = 0.05. Some p are indicated, 
(b): superfluid and vacancy-SS for L = 16 and J/V = 0.1. 



of the states in panels (a) and (b). In fact, in panel (a) 
we performed temperature annealing of the system prior 
to taking the snapshot, while this was not done in panel 

(b) . We find that the defects in (b) reflect the existence 
of a large number of low-energy metastable states, which 
are a direct consequence of the long-range nature of the 
interactions, and will be of relevance for experiments. 

Supersolid and stripe phases are shown in panels 

(c) and (d), respectively. In panel (c), pi has been cho- 
sen to realize an extended vacancy-SS region, surrounded 
by a superfluid. We notice that here a finite T = 0.1J 
has been chosen, compatible with the existence of the 
SS phase (see below). The density-distribution in the 
vacancy-SS looks similar to the ordered CB phase, even 
without annealing. Self-annealing is in fact here enabled 
by the (small) superfluid component of the SS phase. 
Small coherence peaks will be present in time-of-flight 
experiments, allowing for a clear determination of this 
phase. Panel (d) shows a disordered stripe-phase at 
the center, surrounded by an extended Mott-shcll with 
p = 1/4. The disorder in this case is a result of both 
finite T/J = 0.25 and the fact that the stripe solid is 
less robust towards quantum and classical fluctuations 
compared to the CB and star ones. 

These exact results for 51 ^ confirm that the 
phase-diagram Fig. [1] is key to predict and interpret ex- 
perimental observablcs, assuming a local density approx- 
imation. In the remainder of this work, we provide more 
details on the phases described above. 

Incompressible phases: For each Mott lobe 
in Fig. []Ja), the solid order is characterized by 
a finite value of the structure factor S{k) = 
J2 r r ' CX P [*k ( r — r ')] {n r n r i) /N, with k the reciprocal 
lattice vector for each solid. For the CB, stripe, star 
solids, this is (7r,7r), (±27r/3, 2ir/3), and (n,0) [or (0,7r)], 
respectively. The boundaries of the Mott lobes have been 
calculated from the zero momentum Green function (see 
e.g. [H[), for linear system sizes up to L = 20 (CB and 
star lobes), and from p(/x)-curvcs with sizes up to L = 24 
(stripe lobe). 

Interestingly, we find evidence for incompressible 
phases in addition to those corresponding to the lobes 
in Fig. [TJ This is shown in Fig. |3(a,b) , where the particle 



density p is plotted as a function of the chemical potential 
p, for J/V = 0.05, and 0.1, respectively. In the figure, a 
continuous increase of p as a function of p signals a com- 
pressible phase, while a solid phase is characterized by a 
constant p. The main plateaux in panel (a) correspond 
to the Mott lobes of Fig. [TJ while the other steps are 
incompressible phases, with a fixed, integer, number of 
particles. This progression of steps is an indication of a 
devil's-like staircase in the density, the latter being fully 
realized in the classical limit of zero hopping. 

Since the simulations are necessarily performed for 
finite L (with periodic boundary conditions) and T, 
only the lobes with comparatively short periodicity (e.g., 
p =1/2 and 1/4) and sizeable gaps will be resolved in the 
calculations. Consistently, we find that determining the 
groundstatc configuration for each DS-step directly from 
the simulation is often challenging, since: i) for many 
rational fillings [e.g. p = 7/24 in panel (a)] it would re- 
quire to consider system sizes (much) larger than those 
considered here, and ii) the long-range interactions de- 
termine the presence of numerous low-energy metastable 
states [IBj], which for finite T can result in the presence 
of defects or in disordered structures. However, we note 
that the practical relevance of possible Mott lobes with 
large periodicity in the DS region is somewhat limited, 
since they will most likely not be observable in experi- 
ments with trapped molecules, as shown above. 

Supersolid phases: A very different situation is 
shown by the behavior of p(p) immediately above the 
star plateau of panel (a), and below the CB plateau of 
panel (b) in Fig. [3j Here the density grows smoothly 
with increasing p, signaling a fluid phase, and we find 
that the superfluid stiffness p s = T(W 2 ) is finite, with 
W the winding number. Remarkably, and in contrast to 
previous studies with shorter-range interactions, we find 
that the appropriate structure factor 5(k) corresponding 
to these solids remains finite for chemical potentials pi be- 
low and above each Mott lobe, signaling the existence of 
both vacancy-induced and particle-induced supersolidity, 
with no indication of phase-separation. We measured su- 
persolid behavior around both Mott lobes with p = 1/2 
and 1/4, for J/V > 0.05 and 0.067, respectively. As 
an example, Fig. Ufa) shows the case of vacancy-induced 
supersolidity for J/V = 0.2. For an extended range of 
densities, both the superfluid stiffness p s and the static 
structure factor S'(7r,7r), are finite and size- independent. 
We find evidence that the SS melts into a superfluid via a 
second-order Ising-type quantum phase-transition. This 
is shown by the crossing of the curves in panel (b) of 
Fig.gJ where we plot S(k)L 2 ^" as a function of density 
(the critical exponents 2/3 /v = 1.0366(8) correspond to 
the three-dimensional Ising universality class [l(|). The 
supersolid behavior persists for smaller J/V^-ratios, how- 
ever p s tends to decrease with J/V. While the results 
above point to a generic mechanism for solid/liquid tran- 
sitions in 2D with the SS as the intermediate phase fl7j |. 
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FIG. 4: (color online) Vacancy supersolid for J/V = 0.2: 
(a) p s (empty symbols) and S{tt,h) (full symbols) vs. p, for 
L — 8, 12, 16 and 20 (diamonds, squares, dots, and triangles, 
respectively); (b) S(k)L 2B/v vs. p, with 2/3 /V = 1.0366. The 
crossing indicates an Ising-type second-order transition. 



we have not found evidence of SS phases below and above 
the stripe lobe. This may be due to an extremely low T 
for superfluidity here (we have used T m ; n = J/3L), or to 
the strong ID character of the stripe lobe. 

Finite-T: We studied the melting of the SS into a 
normal phase with increasing T, for the case of vacancy 
supersolidity below the CB solid, with J/V = 0.1. Fig- 
ure [5] shows p s and S(n, n) vs T. The melting of the 
SS proceeds through two successive transitions. First, 
SS melts into a liquid-like phase reminiscent of a liq- 
uid crystal, with zero p s and finite S(Tr,n). The drop 
of p s for T ~ 0.1J in Fig. [5] signals a transition of 
the Kosterlitz-Thouless type, with critical temperature 
Tkt = 7rp s fi. 2 p/2m, and m = 1/2 Jo 2 . Upon further in- 
creasing temperature, S(jr, it) drops to zero for T ~ 2.6 J. 
In panel (b) we show that this is consistent with an 
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FIG. 5: (color online) Finite-T melting of SS at J/V = 0.1: 
(a) p s (empty symbols) and S(ir, n) (full symbols) vs. p, for 
L — 8, 12, 16 and 20 (diamonds, squares, dots, and triangles, 
respectively); (b) S(k)L 20/u vs. p, with 2/3/v = 1/4. 



Ising-type transition, by plotting the expected scaling for 
S(tt,tt) in two dimensions (here, 2/3/v = 1/4). 

In conclusion, we have shown that polar molecules 
on a square lattice will realize exotic solid and supersolid 
quantum phases under realistic experimental conditions. 

Note added: While completing the present work, we 
became aware of a simultaneous, independent study on 
a triangular lattice with f2 = 0, see Ref. [Tsj . 
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